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ABSTRACT 

The  relation  of  the  harmonic  spectral  densities  Pm(k,uj) 

of  wall  pressure  on  a  cylinder  in  incompressible  parallel 

flow  to  pure  turbulent-stress  sources  at  highly  nonconvective 

(low)  wavenumbers  ( j  k|  Uoo/uj<<1)  is  formulated  for  arbitrarily 

large  values  of  |k|a  and  j  k | 6 ,  as  needed  in  application  to 

towed-array  noise.  For  the  axisymmetric  component  of  actual 

concern,  the  result  to  zero  order  in  kU  /a)  is  used  in  con- 

junction  with  model  properties  of  the  source  spectra  assumed 

in  earlier  work  to  characterize  the  dependence  of  the  wall- 

pressure  spectrum  p°  (k,<o)  where  |k|6>l.  Where  |k|a<l,  the 

zero-wavenumber  result  previously  found  for  |k|6<l  still 

applies,  and  where  |k|a>l,  an  additional  factor  ~(ka)2 

enters.  The  results  tend  to  confirm,  within  the  uncertainty 

of  the  source  modeling,  the  kind  of  multiplicative  dependence 

of  p  (k,  oi)  on  ka  embodied  in  the  current  semiempirical 
0 

model. 
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1.  INTRODUCTION 

Theoretical  modeling  of  the  spectral  density  p  (k,<jj)  of 

o 

the  axisymmetric  component  of  wall  pressure  on  a  cylinder  in 

parallel  flow  was  pursued  in  Ref.  1.  This  was  done  in  a 

manner  that,  though  still  essentially  dependent  on  ad  hoc 

modeling  of  certain  spectra  of  the  turbulence  field,  provides 

a  more  uniform  framework  and  guide  than  that  used,  for 

example,  in  Ref.  2.  Specifically,  the  subject  spectrum  P  (k,u>) 

o 

was  related  to  spectra  of  turbulent  velocity  products  (or 
Reynolds  stresses) ,  rather  than  partly  to  these  and  partly 
to  a  spectrum  of  the  wall  normal  component  of  velocity  (first 
power,  not  a  velocity  product).  This  new  relation,  in  general, 
exhibits  rather  than  obscures  the  wave  structure  inherent  in 
the  pertinent  dynamic  equations.* 

The  explicit  relation  of  wall  pressure  to  turbulent  stress 
spectra  applied  in  Ref.  1  was  derived  in  Ref.  3  for  a  certain 
domain  of  low  wavenumbers,  defined  roughly  by  the  conditions 

|k|Uoo/ux<l,  |  k  |  6<  <  1  ,  (1-1) 

where  6  is  a  boundary-layer  thickness  outside  which  fluctuations 
are  considered  negligible.  The  main  purpose  of  the  present  memo 


*At  the  same  time,  it  is  recognized,  the  wave  structure  has  no 
direct  relevance  at  subconvective  wavenumbers,  and  the  convec¬ 
tive  domain,  on  the  other  hand,  cannot  be  dealt  with  explicitly 
without  invocation  of  a  closure  assumption  in  some  form. 
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is  merely  to  extend  the  relation  to  a  domain  omitting  the 
second  of  conditions  (1)  and  to  obtain  generalized  results  for 
this  extended  domain  by  use  of  the  same  type  of  source  model 
principally  entertained  in  Ref.  1.  In  much  of  the  parameter 
domain  of  interest  with  regard  to  towed-array  noise,  one  has 
|k|Uco/u)<<l  (i.e.  subconvective  wavenumbers)  but  |  k |  6> >  1 , 
so  that  the  present  extension  is  needed. 

In  Ref.  1,  notwithstanding  the  first  of  restrictions  (1), 
the  relation  for  wall  pressure  was  applied  with  source  spectra 
modeled  as  though  results  should  be  applicable  even  in  the 
convective  domain.  This  was  done  to  furnish  possibly  indicative 
comparison  with  extant  s an i empirical  models  of  the  wall  pressure 
that  encompass  the  convective  domain.*  Source  spectra  are 
likewise  modeled  here  in  this  same  generalized  form. 

In  Sec.  2  the  perturbation  expansion  of  wall  pressure  based 
on  only  the  first  of  conditions  (1)  is  formulated,  and  the 
zero-order  result  given.  In  Sec.  3  results  for  the  wall-pressure 
spectrum  P  (k,<jo)  are  obtained  in  this  zero-order  approximation 
from  the  intended  type  of  source  model,  and  the  implication 
drawn  for  current  modeling  of  towed-array  wall  pressure.  In 
an  appendix  the  next-order  contribution  to  wall  pressure  is 
derived. 


*In  the  "mixed-source"  formulation  of  Ref.  3  (based  on  the  Poisson 
equation  for  pressure  in  the  Lighthill  formulation  for  incom- 
pressive  flow) ,  the  formal  relation  of  pressure  to  sources  is 
readily  available  without  even  the  first  of  restrictions  (1) , 
and  the  problem  of  treating  the  critical  layer  via  closure 
assumptions  that  enters  the  "pure-source”  formulation  when 
applied  to  the  convective  domain  is  absent.  Presumably,  however, 
this  absence  is  no  virtue  but  implies  sweeping  the  dynamic  structure 
totally  under  the  rug  (cf.  preceding  footnote). 
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2.  THE  RELATION  OF  WALL  PRESSURE  TO  FLUCTUATING  STRESSES  AT 
WAVENUMBERS  WELL  BELOW  CONVECTIVE  BUT  POSSIBLY  ABOVE  THE 
RECIPROCAL  DIAMETER  OR  LAYER  THICKNESS 

In  Ref.  3  a  formal  solution  was  obtained  for  the  wall 
pressure  in  inviscid  incompressible  flow  along  a  cylinder.* 

The  subsequent  expansion  procedure  was  valid  provided  that 
conditions  (1-1)  are  satisfied.  Beginning  with  the  same  formal 
solution  for  the  wall  pressure  (Ref.  3,  Sec.  5),  a  modification 
of  the  scaling  used  in  that  development  can  be  made  which 
removes  the  second  restriction  on  k.  The  resulting  expansion 
in  powers  of  e  =  kU^/ui  of  the  transformed  pressure  amplitude 
Pm(k,w) ,  though  somewhat  more  cumbersome,  is  valid  for 
arbitrary  values  of  k5 .  Since  the  method  to  be  presented 
closely  parallels  the  procedure  of  Ref.  3,  only  a  brief  outline 
of  the  results  will  be  recorded  here.  Furthermore,  explicit 
solutions  will  be  given  only  for  the  axisymmetric  component 
(m=0)  of  the  wall  pressure,  though  the  method  is  applicable 
for  any  circumferential  mode. 

The  transformed  pressure  amplitudes  p  (r,k,w)  of  the  cir- 

cumferential  modes  e  are  related  to  source  function  §  , 

m 

which  contains  components  of  the  fluctuating  stress  tensor,  via 
the  ordinary  differential  equation 


*This  was  done  also  for  slightly  compressible  flow,  but  that 
aspect  is  not  of  concern  here. 
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The  explicit  form  of  is  given  in  Ref.  3,  Eqs.  37  and  38. 

The  formal  solution  of  (2-1) ,  subject  to  the  boundary  conditions 

l£m(r)  |  <  °°  at  r  =  « 

and 

I  ft '  (r)  |  =  0  at  r  =  a  , 

1  Mn 

the  cylinder  radius,  may  be  obtained  in  a  straightforward  manner. 
When  evaluated  at  r=a,  the  axisymmetric  component  becomes 

00 

Po(a,k^)/p  =  A“1rc(a)|  dr  (Trr-Txx)  k2r<f>Ec  +  (Trr+T0e>  <*> 'Ec 

-  <Trr-Teeti2teTrx»*E;(r)]  <2'2) 

where 

r  (s)  a  <f>F (r)/u 

Ec(r)  a  4>e  (r)  /go 

A  *  r<Ecrc-rcEc> 

4>  a  (1-kU/oj)”1  , 

and  E  and  T  are  the  solutions  of  the  homogeneous  form  of  (2-1) , 
satisfying  the  boundary  conditions 

|E  (r)  |  <  00  at  r  =  00 

|  T*  (r) |  =  0  at  r  =  a  .  (2-3) 
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A  low-wavenumber  expansion  for  p  requires  corresponding 

o 

expansions  for  the  homogeneous  solutions  E  (r)  and  r  (r) . 

c  c 

These  may  be  obtained  more  conveniently  by  considering  a 
self-adjoint  form  of  (2-1) , 


ka+ai 

r2  r%j> 


g_  =  0 
m 


(2-4) 


where 


g  (r)  *  p  (r)rv(u)-kO)  . 
m  m 


Adopting  the  scales  £  =  r/6,  U  =  U^VfC),  e  =  kU^/ui  and  s  =  k6, 
(2-4)  becomes 


LZa .  fs^sa-v<  + 

d£2  t  £2 


£V*  .  £V"  .  2ezV 2  _  _  „ 

+  i=w  +  r*.  s*  -  o  ■  (2-5) 


(1-eV)' 


which,  when  expanded  to  second  order  in  e,  yields 

— fs.  -  |s2+2-ZiZi  +  e(v"+~)  +  e 2 [ 2V '  2+V lV"+~)  ]  j  g  =  0. 

dC2  52  g  *  I 

(2-6) 

Inserting  a  power  series  representation  for  gm, 

gm  -  3.  +  +  c‘g,  +  ••• 

results  in  the  following  hierachy  of  equations 

e°:  g"  -  [ s2+ (m2-l/4) /C2 ]g  =  0 

0  0 

£*:  g"  -  [s2  +  (m2 — 1/4 ) /5 2 ]g  =  g  (V"+V'/S> 

1  10 

e2"  g"  -  [s2+(m2-l/4)/C2lg  =g  (V”+V'/Q+g  [2V '  2+V  (V"+V  */£)  ] 
2  2  10 

(2-7) 
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The  boundary  conditions  identifying  the  two  solutions  E  and 
F  must  also  be  expanded  in  powers  of  e.  The  wall  condition 
p^(a)  =  0,  corresponding  to  7,  can  be  written 

[U-eVUnr^U)]  =  0  at  C  =  a/6  e  Cq  , 

which,  when  expanded,  supplies  a  boundary  condition  for  each  of 
the  equations  in  (2-7) 

e°:  g'  (Co)  =  (1/2). C_1g 
0  0  0 

e1:  g1  (C  )  =  g  V’  (C  )  +  k_1g  (2-8) 

1  0  0  0  2  0  1 

e2:  g'(C  )  =  gV'(C  )  +  k-1g  etc. 

20  1  0  2  0  2 

The  boundedness  condition  on  p  (r)  at  large  distances  from 

m 

the  cylinder  translates  immediately  into  the  condition 

-L 

C  9m^<0°  as  Ifc  more  convenient,  however,  to  form 

an  equivalent  condition  which  is  applied  at  £  =  1.  Outside  the 
boundary  layer,  where  £>1,  (2-5)  simplifies  to 

d2g 

- £  -  [ s2+ (m2-l/4)/Cz]g  =  0  , 

dC2  m 

which  has  as  its  finite  solution 

gmU>  =  ac^^ ( |  s|  C)  .  (2-9) 

The  solution  for  C<1  must  merge  with  (2-9)  in  such  a  way  that 
g^ ( C ) /gm ( C)  is  continuous  at  C=l,  or,  equivalently. 
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dg ( 5 ] 

m 

d? 


K  (  i 

m 


S1?)  =  *m( 


5,+aK( 


?)  ]  at  4=1  • 


Inserting  the  expansion  for  gm  results  in  a  sequence  of  boundary 
conditions  to  be  applied  at  ?=1, 


dg.(C) 

d? 


( I  s  j  ?) 


gi(C)  ^IKm(  I  S)+^Km(  I  sl 


(2-10) 


on  the  e1  level  of  the  hierarchy  given  in  (2-7) .  This  method  of 
selecting  the  bounded  solution  corresponding  to  E  proves  to  be 
more  tractable  when  the  higher-order  (inhomogeneous)  equations 
of  (2-7)  are  considered. 

The  zero-order  ( e 0 )  solutions  for  the  axisymmetric  mode 
(m=0)  may  be  readily  obtained  by  solving  the  first  equation  of 
(2-7)  subject  to  (2-8)  and  (2-10)  .  The  solution  bounded  at  °°  is 


g  (?) 

0 


B  (  |  s | ?)  , 

0  0 


(2-11) 


and 


g  (?)  =  A  ?' 

0  0 


I  ( 

0 


I'(  s  ?  ) 

s|  ?) — 2 - 2- 

K  *  ( J  s  1 C  ) 

0  0 


Ko( 


S  I  ?) 


[< 


satisfies  the  proper  wall  condition  at  ?=?  .  To 

0 

with  the  normalization  A  =B  =1,  the  functions  E„ 

00  c 

become 


) 1  =  d/d?]  (2-12) 

this  order,  and 

and  T  in  (2-2) 
c 


E  (r)  *  K  ( |k|r) 

C  0  I ’ ( | k| a) 

:  (r)  =  I  (|k|r) — 2 - K  ((]4r)  .  [(  )'  =  d/dr" 

0  K'(ikja)  0 
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This  implies  the  following  zero-order  (ins)  approximation  for  the 
axisymmetric  component  of  the  wall  pressure 


p  (a,k,u i)/o 
0 


1 

j  ka | K  ( | k i a) 


[  dr[  (f  _f  )k2rK  (|k 
j  rr  xx  o 

a 


r) 


+ 


(Trr+Tee) (k/^)u’ (r) Ko (Ikir)  + 


('?rr'^ee+2ikr^rx)  f  k  i  K x  C I  k |  r )  ]  .  (2-13) 

In  the  limit  k5  (and  ka)-*-0,  Eq.  2-13  attains  the  proper  form 
(Ea.  58,  Ref.  3).  It  provides  an  expression  for  the  wall  pressure 
(to  zeroth  order  in  )  having  an  extended  domain  of  validity. 

In  the  appendix  it  is  indicated  how  an  additional  order  of 
accuracy  (in  z)  may  be  achieved.  The  following  modeling,  however, 
will  employ  only  the  zero-order  result  (2-13) . 


3.  THE  WALL-PRESSURE  SPECTRUM  BASED  ON  ASSUMED  PROPERTIES 
OF  SOURCE  SPECTRA 

We  wish  to  use  the  relation  (2-13)  of  wall-pressure  amplitude 
in  conjunction  with  specific  assumed  properties  of  the  cross-spectra 
Ef  jk^  ( r » r  ' »  k , u))  of  fluctuating  stress  components,  corresponding  to 
the  product  TV*  (r ,  k,  us)  T^  (r  ' ,  k  ’  ,w),  to  infer  a  resulting 

expression  for  the  wall -pressure  spectrum  P  (k,a))  and  examine 

0 

its  limiting  dependences.  Specifically,  we  assume  a  model  of 
the  type  discussed  in  Ref.  1  in  connection  with  Eqs.  16  to  19, 
and  25.  In  this  context,  the  model  is  exemplified  by  the  form 
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Eijk£  (r,r  '  ,k,ui)  =  cij£mv*  (a‘/rr  1  )L(r)L  (r  ' )  (l+ak+r)  (l+ak+r '  ) 

x  S* (k+v) S (k+y ' )  .  (3-1) 

(r=a+y,r '=a+y 1 ) ,  where  the  scale  length  L  is  taken  to  be  given 
by  the  mixing-length  model  form  (Eg.  25) 

u 

L(r)  ~  (ar)  £n(r/a)  ;  (3-2) 

k+  =  [  (u>-uck)  2/ (hv*)  2+k2 ]**  ,  (3-3) 


ct~l,  u^U^,  the  cj_j£m  ar©  dimensionless  constants  (~1),  and 
S(k+y)  is  a  well-behaved  function  that  is  '1  for  k+y<l  and 
decreases  rapidly  for  k+y>l.  Though  not  expressed  in  the 
notation,  S(k+y)  is  assumed  negligible  for  y>5  independently  of 
k+  and  likewise  to  vanish  appropriately  within  the  viscous  sub¬ 
layer  as  yv*/v-*-0 . 


P  o  (k,  cu)  inferred  by  reference  to  (2-13),  as  usual  (see  Ref.  1), 

is  expressible  as  a  sum  of  a  number  of  double  integrals  over  r 

and  r'.  Under  assumption  (3-1),  each  of  these  is  simply  a  square 

of  a  single  integral.  To  simplify  expressions,  temporarily 

imagine,  for  example,  that  c. ..  in  (3-1)  also  had  the  factorable 
*  ijx-m 

form  cj_jck£*  In  (2-13)  the  mean  velocity  gradient  outside  the 
sublayer  may  be  approximated  as  in  Ref.  1,  Eg.  24,  by 


U'(r)  =  Av*/r£nr  (A=2.5),  r-a>6v/v*.  (3-4) 

The  spectral  density  obtained  from  (2-13),  with  tentative 
neglect  of  source  strength  within  the  sublayer,  is  then  given  by 

P,  (k»u,)  »  c  £v*a 2 1 2 ,  (3-5) 
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TTa'K-TikTa)  j  dyS(k+y)[l+ak+a(l+y/a)]- 
1  6v/v* 


(1+y/a)  *£n  (1+y/a)  {  (c  -c  )  k:  a  (1+y/a)  K  'jk|a(l+y/a] 

J-  L  XX  q 


+[crr-cQ6  +  i2ka (1+y/a) crx3 [ k [ K ^  ( [ k | a (1+y/a) ) 


+  ^crr+cee ^ k (Av*/wa)  (1+y/a) -1[2n  (1+y/a) ]_1Kr  ( jk|a(l+y/a) ) } 

( 3-5a) 

The  essential  aspects  of  the  dependence  of  (3-5a)  and  hence 

P  (k, w)  under  present  assumptions  are  demonstrated  by  con- 

o 

sidering  the  opposite  limits  ka<_<l  and  ka>J>l  where  in  both 
instances  ko>^>l.  [The  case  where  ka<<l  and  k5<<l  was  already 
dealt  with  in  Ref.  1  and  corresponded  to  approximating  {  ;  in 
( 3-5a)  by  the  dominant  limiting  terms  (Crr-Cpe )  j  k  |  K  j  (  j  k  j  r ) -*■ 

(crrc00)r“1*3 

First,  when  |k|a<<l  and  ! k | 6 > > 1 ,  since  k+>|k|  and  the  factor 
S(k+y)  cuts  off  contributions  to  the  integral  at  y-k"1,  one  has 
for  values  of  y  that  may  contribute  appreciably: 

| k i (a+y) < | k j  (a+k" 1 ) <  k  ( a+k" 1 ) -1 .  It  is  thus  not  a  bad  approxi¬ 
mation  to  regard  |k|  (a+y)  in  K ^ ( i k j  (a+y))  as  small.  Hence, roughly 
the  same  result  as  the  earlier  one  for  ] k | 6< < 1  is  again  obtained.* 


*That  result  corresponds  to  replacing  the  lower  limit  6v/v*  in 
(3-5a)  by  zero  and  hence  loses  validity  where  6k  v/v*>l.  The 
same  is  true  of  the  result  given  for  | k i a > > 1 . 
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Next,  when  ,k,a>>l,  the  K^(,k. (a+y) )  may  be  approximated 
for  the  entire  range  of  integration  by  their  limiting  forms 
for  large  real  argument.  Eq.  3-5a  then  becomes 

6 

I  -*■  a  dyS  (k+y)  exp  (-j  k|  y)  [l+ak+a  (1+y/a)  l”*5 

6v/v* 


*  { in  (1+y/a)  [  I M  a (crr+i2crx~cxX) -(1+y/a)  "J  (cQe-crr)  ] 


+A(v*/aia)  (1+y/a)  "2  (crr+cQe)  } 


(3-6) 


for  k<0,  respectively.  Here  Av*/coa<<l  (even  if  k-Lo/U^,  since  ]k|a>>l).* 

Barring  unexpected  cancellation  among  c.  .'s,  the  term  11  c  +c,,as 

ij  rr  vc 

well  as  that  «  cee~crr  may  ne9lected*  Further,  S(k+y)  cuts 
off  contributions  (since  k  >>5_1)  above  y~k" l< I k I ” 1 ;  for 
contributing  y,  then,  y/a<l/ | k j a< <1 .  Expression  (3-6)  then  reduces, 
apart  from  a  factor  of  the  order  of  unity,  to 


I  ~  <crr+i2crx-cxx)|kik;5/2  ~  I k“5/2 


and  (3-5),  in  turn,  to 


P  (k,w)  ~  o2v2 (ka) 2k”5 . 
0  + 


(3-7) 


Relaxing  the  assumption  that  the  c^..^  in  (3-1)  have  a  product 
form  and,  in  fact,  generalizing  S* (k+y) S (k+y ' )  to  become 
S* j (k+y) Sk£ (k+y ' )  (with  similar  properties)  will  evidently  have  no 


*More  generally,  for  j  k | a  >  > 1 ,  Eq.  2-13  becomes 

p  (a,k  ,w)/c  -  a-1  dy  (1+y/a)  _i5exp(  - :  k  y)  [  (c  +i2c  -c  )  k’(a+y) 
c  rr*  xx 

*  r- 

+c  -cC£.+  lc_  +cai,)U'  (a+y)/-] 

iT  i  v-  v  Li  -  - 


11 


Tech  Memo  23 


Chase  Inc . 


effect  on  the  functional  dependences.  Likewise,  these  dependences 
would  prevail  for  a  class  of  spectral  functions  somewhat  more 
general  than  the  product  form  of  Eq.  3-1. 

In  the  limit  of  small  |k|  (a+6)  ,  p  (k,oj),  it  is  noted,  is 

o 

dominated  by  the  contribution  from  the  source  spectral  density 

related  to  Trr-Tgg,  and  in  the  large-|k[a  limit  to  that  related 

rather  to  T  +i2T  -T 

rr—  rx  xx 

The  form  (3-7)  agrees  with  the  form  assumed  by  the  current 
semiempirical  model  for  P  (k,u>)  (e.g.,  see  Ref.  1,  Eq.  9a) 

in  the  pertinent  limit.  The  present  results  for  functional 
dependence  on  ka  result  from  the  form  of  the  homogeneous 
operator  on  pressure  in  the  underlying  equation  and  from  the 
cylindrical  geometry.  As  it  turns  out  (and  as  one  might  hope) , 
these  results  based  on  an  underlying  equation  consistently 
involving  velocity-product  sources  yield  no  departure  from  those 
obtained  in  Ref.  9  based  on  a  Lighthill  formulation  involving  mixed 
sources.*  Hence,  theoretical  modeling  still  does  not  suggest 
that  a  trial  form  different  from  our  current  wall-pressure  model 
given  in  Ref.  1,  Eq.  9a  (with  term  linear  in  ka  included)  is 
likely  preferable. 

The  question  of  possible  entrance  of  another  length  scale, 
as  discussed  in  Ref.  1,  and  of  the  numerical  values  of  parameters 
in  the  model  have  not  been  further  addressed  here. 


♦Since  Ref.  2  was  written,  the  assumed  form  of  the  detailed 
dependence  on  k+  has  evolved  to  a  somewhat  different  form, 
as  discussed  in+Ref.  1,  Sec.  3. 
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APPENDIX.  EXPANSION  OF  WALL  PRESSURE  TO  NEXT  ORDER  IN  kU  /u 

00' 

WITH  k5  AND  ka  ARBITRARY 

The  components  necessary  for  an  0(C)  correction  to  p  ,  given 

o 

by  (2-2),  are  calculated  in  this  appendix.  This  requires  the 
solution  of 

g" ( C ) - ( s2-l/4£ 2 ) g  =  g  [v ” ( S) +V 1 /C]  (A-l) 

1  10 

subject  to  the  boundary  conditions 

g'(5  )-kr'g  (C  )  =  g  (C  )V' (5  )  (A-2) 

l  o  *  o  l  c  oo  o 

and 

g '  ( C )  K  ( |  s  j  C)  -  g  ( C)  C  4k  (|s)C)+^k  (|s|£)]  at  5  -  1  , 

io  i*o  d  o 

(A-3) 

which  yield  Fc  and  Ec,  respectively. 

The  general  solution  to  (A-l)  may  be  written 

r 

g  (C)  =  A  f  +B  j-  -  [p  (C)^  ( t)  — ij;  (C)'i'  (t)]g  (t)  [v"  (t)+V'/t]dt 

1  1  1  1  2  J  2  1  1  2  0 

„  5.  „ 

where  <p  =■  (|s|C),  y  =  £  K  (|s|£)  are  the  independent 

10  2  0 

solutions  of  the  homogeneous  form  of  (A-l).  The  ratio  of  the 

constants  A  ,  B  is  determined  by  a  normalization  condition,  and 
i  i 

that  ratio  is  then  fixed  by  the  application  of  either  (A-2)  or 
(A-3) . 

Consider  first  the  boundary  condition  (A-2)  ,  corresponding 

to  the  T  solution.  With  g  (t)  given  by  (2-12),  ( A— 4 )  may  be 

0 

written 
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C_i5g  U)  =  I  ( |  s| 5) 

1  3 


,  I '  (  ;  s  i  l  ) 

A  +  j  K  (  j  s  i  t)  ( I  (  |  s  |  t ) 2 - 2 —  K  (jsjt)]  (tV'(t))'dt 


i  J 
4 


k  1  ( ;  s  j  c  ) 
0  0 


+  K  ( |  s| c) 

o 


B  -f  I  (  |  S  j  t )  / 1  (  [  S  ]  t) - S - 

1  l  0  l  0  K  '  (  j  s  j  5  )  ° 

c,  \  oo 

o 


I'(|s|5  )  v 

-° - S—”  (  |  S  |  t)  J 


(tv  (t))  'dt 


Applying  the  boundary  condition  (A-2)  yields  the  following  relation 


between  A  and  B  , 

1  1 


A  I '  (  !  s  I  4  )+BK'(|s|£  )  *  -V 1  ( C  )/C  K  '  {  |  s  I  c  ) 
10  010  0  000  0 


The  normalization  condition  g  (5  )  =  0,  a  choice  for  convenience 


l  o 


only,  provides  the  additional  condition 


0  «  A  I  (|s|  ?  )  +  B  K  (|  s|  5  )  . 

i  o  o  ioo 

With  A  ,  B  so  determined,  the  first-order  correction  g  (5)  becomes 
ii  i 

,  -V(5  ) 

4  \  (4)  = - - -  [K  (|s|4  )I  (  I  S  ]  4 )  “I  (  ,  sl  4  )K  (|s|5)] 

1  K*(|S|C  )  ° 


-j  (tV'(t))' 


K  <  |  s|  C)I  (Isjt)-K.  ( I  s  1 1)  I  (|sU) 

0  0  oc 


l’(|s|£  ) 

x  ( !  S 1 1)  — 2 - L  K  ( |  s  i  t)  dt 

0  K  '  ( !  s| 4  )  0 

oo  - 


or, equivalently ,  the  0(5)  correction  to  rc(r)  in  (2-2),  denoted 

by  T  (r) ,  is 
i 

r  ( r )  =  Z.Vf-v:  ■  '  [k  ( i  k  |  a )  I  ( |  k  i  r )  -K  (  I  k  |  r )  I  (  !  k  |  a )  ] 
i  JvViKiaio  o  o  o 


frd 

3T(tv’ 

'a 


! 

I’(lkia) 

<t) ) 

K_  (  1  k ; r )  I,  ( 1  k,  t)-sr  ( '  k!  tEr  (|kir) 

I  ( !  k  t)  — - - k  .  ( 1  k;  t) 

:  K'(ikia)  ; 

c 

( A— 5 ) 
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where  V(r)  =  Ulr)/!^.  It  is  readily  established  that  in  the 
limit  kc  (and  ka)-*-0,  r  (r)  given  by  (A-5)  tends  to  v(r) 
so  that  to  0(e) 


rc  (r)  *  l+eV(r) , 

in  agreement  with  the  results  of  Ref.  3. 

In  a  similar  manner,  the  first-order  correction  to  Ec(r) 
is  calculated  by  subjecting  the  general  form  (A-4)  to  the 
appropriate  boundary  condition  (A-3)  and  using  <3g(i)  as  given 
by  (2-11).  The  normalization  condition  =  0  again 

provides  the  relation 


0  =  A  I  (|s|5  )  +  BK  (|S  C  )  , 
10  0  10  0 


and  (A-3)  in  turn  yields 


A  =  -  K2  ( |  s|  t)  (tV'  (t)  )  'dt. 
i  J  o 


The  appropriate  solution  is 


(O  =  -IQ  (|sU) 


K2  ( 1  s  1 1 )  ( tv  ’  ( t) )  '  dt 


LI 


-K  ( I  s|5  ) 
0  0 


1 


•  I  (is|C  , 

I  o  ( I  s  i  t)K  0  ( I  s  1 1)  (WM’dt-^j^l  K20(|s|t)(tV')'dt 

0  Or 

n  ^  _ 


which,  when  written  as  a  function  of  r,  denoted  by  E ^  (r) , 
provides  the  first-order  correction  to  E  (r) , 


E_ (r )  =  K  (kr)  +  ZE  (r)  , 

c  c  1 
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with 


E  (r) 

i 


-I  ( j  k  |  r )  {  k2  ( |  k  1 1)  (tv')  'dt 


V  ,  I.(|k|a)' 

5 

f  „  , 

K  ( i  k | r) 

0 

[  K 2  ( |  k|  t)  (tV')  'dt 

L  a  u  a  J 

(A-6) 

The  limiting  form  of  (A-6)  when  k6  -*•  0  properly  reduces  to  the 
formula  given  in  Reference  3,  where  Ec(r)  was  calculated  with 
kf  assumed  small  at  the  outset. 
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